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ABSTRACT 

The Magellanic Stream and several high velocity clouds have now been detected in optical line emission. The 
observed emission measures and kinematics are most plausibly explained by photoionization due to hot, young 
stars in the Galactic disk. The highly favorable orientation of the Stream allows an unambiguous determination 
of the fraction of ionizing photons / esc which escape the disk. We have modelled the production and transport 
of ionizing photons through an opaque interstellar medium. Normalization to the Stream detections requires 
/esc ~ 6%, in reasonable agreement with the flux required to ionize the Reynolds layer. Neither shock heating nor 
emission within a hot Galactic corona can be important in producing the observed Ha emission. If such a large 
escape fraction is typical of L» galaxies, star-forming systems dominate the extragalactic ionizing background. 
Within the context of this model, both the three-dimensional orientation of the Stream and the distances to high- 
velocity clouds can be determined by sensitive Ha observations. 

Subject headings: Galaxy: halo - galaxies: individual (Magellanic Stream) - galaxies: ISM, intergalactic medium 
- cosmology: diffuse radiation - instrumentation: interferometers 



1. INTRODUCTION 

There is extensive theoretical and observational interest in 
establishing what fraction of the total ionizing luminosity from 
the stellar disk of the Milky Way and other galaxies escapes into 
the halo and the intergalactic medium (e.g., Miller & Cox 1993; 
Dove & Shull 1994; Leitherer et al. 1996). This quantity is also 
of cosmological interest, since if / esc exceeds a few per cent, 
then star-forming galaxies dominate over AGN in producing 
the ionizing background at low redshift (Giallongo, Fontana, 
& Madau 1997). Observationally, it is extremely difficult to 
measure this quantity directly. Leitherer et a/. (1995) used HUT 
observations of four UV-bright galaxies to constrain / esc to typ- 
ically no more than a few percent at the Lyman limit. However, 
Hurwitz, Jelinsky & Dixon (1997) showed that absorption by 
gas within the Milky Way, which was ignored by Leitherer et 
al, raises these limits substantially, to ^ 10%, so that they are 
no longer significant constraints. 

Recently, Weiner & Williams (1996; hereafter W 2 ) have 
detected Ha emission at three points along the Magellanic 
Stream. We demonstrate that this signal is almost certainly due 
to ionizing photons escaping from the disk of the Galaxy. In the 
following section, we present a simple model for the production 
and escape of ionizing photons from the disk. In §3, we use the 
Ha observations to measure the mean ionizing photon opacity 
(f ) of the disk, which determines both the shape and normal- 
ization of the Galactic halo ionizing field. Finally, in §4, we 
discuss the implications for observations of the Stream, high- 
velocity clouds, and the cosmic ionizing background. 



2. GALACTIC PHOTOIONIZATION MODEL 

Galactic disk. We construct an idealized model for the ionizing 
radiation field away from the Galactic plane. The coordinate 
system is shown in Fig. 1. The ionizing stars are assumed to 
be isotropic emitters confined to a thin disk; since we are in- 
terested in emitting regions at distances much greater than the 



stellar disk thickness, the neglect of the finite thickness of the 
disk is reasonable. At coordinates (x o ,0,z o ), a distance R„ from 
an arbitrary patch of the disk dA, the received flux fj (in units 
of erg cm" 2 s" 1 Hz" 1 ) from ionizing disk sources with specific 
intensity £„, through a solid angle dft, is 



fd 



C„ drj = / cos 9 



dA(r ,4> ) 



Rl 



(1) 



where dA = r a dr d<f) and R^ = x 2 + z 2 + r 2 - 2x r cos <j> . The 
normally incident, ionizing photon flux from the disk (in units 
of phot cm" 2 s" 1 ) is 

dA(r ,<j> ) f da(r ,<p )cos9 



<Pd(r ,<l>a)= / n d (r ,<j) ) cost 
J A 



Rl 



Rl 



(2) 

for which «,/ and der are the surface photon density and bright- 
ness, respectively, within each disk element dA. 

Vacca et al. (1996) have compiled a list of 429 O stars 
within 2.5 kpc of the Sun, from which they determine an ion- 
izing photon surface density at the Solar Circle of nd( r 0) ~ 
3.74 x 10 7 phot cm" 2 s" 1 , where r Q ~ 8.0 ±0.5 kpc (Reid 
1993). We assume an axisymmetric exponential disk with 
radial scale length rj, so that rid(r) = n e~ r ' Td . With = 
3.5 kpc (Kent, Dame & Fazio 1991), we thus derive n ks 
3.7 x 10 8 phot cm" 2 s" 1 . The explicit contribution from B stars 
leads to only a small correction. The O star survey region 
includes part of the Local Arm and extends to the neighbor- 
ing Sagittarius-Carina and Perseus spiral arms (see Taylor & 
Cordes 1993). Observations with the Ultraviolet Imaging Tele- 
scope of face-on, mid to late-type spiral galaxies (made avail- 
able by E.P Smith) give us confidence that our extrapolation is 
reasonable. All integrations are performed out to 25 kpc in ra- 
dius {i.e., de Geus et al. 1993; Ferguson et al. 1998), giving an 
integrated photon flux over the disk of 2.6 x 10 53 phot s" 1 . This 
value is consistent with Mezger's (1978) estimate and within 
25% of the rate inferred from COBE/FIRAS observations of 
the Galactic [Nil] 205^m emission (Bennett et al. 1994). 
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Fig. 1. Coordinate system for a point D(x o ,0,z o ) situated above 
an opaque disk, and a point H(Q,y 1 ,z\) within an extended halo. 
The limit of integration r\ for the halo is a function of 6\ . 
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Fig. 2. An illustration of the LMC and the dominant clouds 
in the Magellanic Stream (Mathewson & Ford 1984) projected 
onto the Galactic X-Z plane. The orbit of the Stream lies closer 
to the Great Circle whose longitude is / = 285°. We ignore 
small projection errors resulting from our vantage point at the 
Solar Circle. The angle 5 is measured from the negative X axis 
towards the negative Z axis where 8 = -b (0° < 8 < 90°) and 
<5 = Z?+180° (90° < 8 < 180°). 



For an opaque disk, the patch dA is observed through the in- 
tervening disk interstellar medium (ISM) with optical depth r 
(at a given frequency), such that the observed surface brightness 
da' = e~ T( - r °'^ da. We approximate the disk ISM by a smooth 
distribution. Since the absorbing disk thickness (^100 pc) is 
negligible compared to the height of the Magellanic Stream 
above the Galactic plane, we can assume plane-parallel absorp- 
tion for which T{r ,cf> ) = kR /z , where n is the photon absorp- 
tion perpendicular to the galactic disk. Since for the purposes 
of this paper we are only interested in the total ionizing photon 
flux, and not the detailed spectral shape, we characterize the 
radiation field by the mean ionizing photon opacity f (i.e., the 
value averaged over photon energy). 

Galactic corona. Spitzer (1956) first suggested that the Galaxy 
has a hot corona, with a scale height of several kiloparsecs. Ab- 
sorption line observations of highly ionized species have con- 
firmed the presence of hot gas at kiloparsec heights above the 
Galactic plane (e.g., Savage et cd. 1989; Sembach & Savage 
1992). The presence of an actual large-scale hot corona, as op- 
posed to a patchy distribution of hot gas arising in the disk, is 
more uncertain (e.g., Wang & McCray 1993). To investigate the 
influence of such a corona, we assume a distribution of hot gas 
of the form n(R\) = n c /(\+R\/rl) cm" 3 , where n c is the core 
density and r c the core radius, and a gas temperature T h . 

The halo is assumed to be finite in size with radius R h . The 
normally incident flux on a (plane parallel) cloud face, tangen- 
tial to the radius vector, at a distance R\ from the origin is de- 
scribed by 



//, = f % cos6>i dfl = 2n 



tt/2 



F 6 



Ri Ri ' 



(3) 



for which is the photon emissivity, F is a somewhat compli- 
cated function of R\, 6\, r c , and the (6\ -dependent) radial limit 
of integration r\ , which is given by 



n = 



R„ sin(6»i + sin" 1 [Ri sin6>i /R,]) 



sinful 



(4) 



The photon emissivity £„ was computed with the MAP- 
PINGS 3 ionization code (Sutherland & Dopita 1996) and in- 
corporated bremsstrahlung, line and non-equilibrium cooling. 
The ionizing photon flux on the inner face of a cloud is given to 
10% accuracy (for R , /r c < 1 2) by 



p*(Ri) « 10V 3 r c 



6 



10" 



77+1.3(fli/r c ) 
(\+R\/riy 



1.35 



phot cm s 



where 10 3 n_3 cm 3 is the core number density, R\ and r c are in 
kpc, and £, is the frequency-integrated ionizing photon emissiv- 
ity. For an optically thin Galactic disk, -q = 7r/4; when the disk 
is completely opaque to ionizing photons, 77 = 0, which reduces 
the ionizing photon flux by less than a factor of two. 

Estimates of T,„ r c and n c have come from three independent 
lines of inquiry. Absorption lines towards stars in the Magel- 
lanic Clouds have been used to infer coronal temperatures and 
gas densities; e.g., Songaila (1981) finds T,, ks 0.05-0.3 keV 
(0.6-3.5 x 10 6 K) and n c ~ 3-7 X 10" 4 cm" 3 . Alternatively, 
various groups have attempted to model the Magellanic Stream 
as clumps moving through a hot, viscous fluid (q. v., Moore & 
Davis 1994). An additional constraint is provided by pulsar 
dispersion measures (V m = 25-100 cm -3 pc) along sight lines 
to the LMC/SMC system, and a halo globular cluster, NGC 
5024 (Taylor, Manchester & Lyne 1993); these values imply 
electron densities of order 10 -3 cm -3 . However, the D m col- 
umn arises primarily within the diffuse ionized gas layer in 
the Galactic disk (the Reynolds layer), so that the true elec- 
tron densities must be considerably smaller. Following Moore 
& Davis (1994), we adopt r c = 10 kpc and n_3 = 2, which pre- 
dicts a density of 6 x 10 -5 cm" 3 at R„ = 55 kpc, and assume 
T„ w 0.2 keV « 2.3 x 10 6 K, the virial temperature of the Galac- 
tic halo (for an assumed circular velocity of 220 km s _1 ). This 
density and temperature distribution matches the emission mea- 
sure and X-ray luminosity of the diffuse component determined 
by Wang & McCray (1993). (See also Maloney & Bland- 
Hawthorn 1998.) For this temperature, ~ 0.3-0.4 x 10 -15 
cm -3 s _1 sr _1 , with only a weak dependence on metallicity. 
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arrive at the cloud surface (phot cm -2 s" 1 ), n e is the electron 
density and N/ is the column density of ionized hydrogen and 
helium. 

The emission measure is £ m = J n e ni dl = n e riiL cm" 6 pc 
where L is the thickness of the ionized region. The result- 
ing emission measure for an ionizing flux ip-, is then £ m = 
1.25 x I0~ 2 ip4 cm" 6 pc (= 4.5(^4 mR) where ip t = lO 4 ^- For 
an optically thin cloud in an isotropic radiation field, the solid 
angle from which radiation is received is il = Air, while for one- 
sided illumination, f2 = 2ir. For our disk model, however, J v is 
anisotropic and £1 can be considerably less than 2ir. 



Fig. 3. HI peak intensity map from the HIPASS survey of the 
Magellanic Stream (Putman et al. 1998) centered on two W 2 
fields, where the W 2 beam size and pointings are shown by cir- 
cles. Strong Ha detections fall in regions of high mean column 
density (e.g., MS II A); weak detections (e.g., MS II B) and 
non-detections fall in regions of much lower column density. 
The contours are for brightness temperatures of 0.1, 0.2, 0.4, 
0.6, 1, 1.4, 1.8 K. 

These plasma constraints on the density and temperature im- 
ply that any contribution from a Galactic corona to the ioniz- 
ing photon flux at the Stream is insignificant (eqn. ||). For the 
known velocity of the Magellanic system (ss215 km s ; Lin, 
Jones & Klemola 1995), and for volume-averaged, cloud densi- 
ties of roughly 1 cm" 3 , to produce substantial Ha through shock 
ionization requires halo densities which are orders of magni- 
tude higher than inferred from the pulsar dispersion measures. 

3. PHOTOIONIZATION OF THE MAGELLANIC STREAM 

The Stream lies along a great arc which extends for more 
than 100° (e.g., Mathewson, Cleary & Murray 1974). Fig. 2 il- 
lustrates the relationship of the LMC to the Magellanic Stream 
above the Galactic disk (Mathewson & Ford 1984). The LMC 
lies close to the X-Z plane. We make the assumption that the 
Stream lies along a circular orbit, close to the X-Z plane, origi- 
nating from the Lagrangian point between the LMC and SMC, 
at a Galactocentric distance of 55 kpc. 1 However, most com- 
puted orbits for the LMC-SMC system imply substantial ellip- 
ticity along the Stream (e.g., Lin et al. 1995). 

In order to derive Ha emission measures, we assume an elec- 
tron temperature T f »s 10 4 K, as expected for gas photoionized 
by stellar sources; the Case B hydrogen recombination coeffi- 
cient is a B « 2.6 x 10" 13 (10 4 /7;) a75 cm 3 s" 1 . At these temper- 
atures, collisional ionization processes are negligible. In this 
case, for an optically thick slab, the column recombination rate 
in equilibrium must equal the normally incident ionizing pho- 
ton flux, a B n e Nj = (ft, where </?,- is the rate at which Lyc photons 

'For the LMC, cepheids indicate (m-M) = 18.47 ±0.15, which implies 49.4±3.4 
kpc (Feast 1988). 

2 1 Rayleigh is W 6 /4ir phot crrT 2 s" 1 sr~' or 2.41 X 1(T 7 erg crrT 2 s" 1 sr~' at Ha. 
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Fig. 4. The predicted Ha emission measure along the Stream as 
a function of S (Fig. 2). The vertical axis has units of log(cm" 6 
pc), equivalent to log(Rayleighs) after subtracting 0.48. The top 
curves assume an optically thin Galactic disk with (dashed line) 
and without (solid line) the LMC ionizing field. For the lower 
curves, the solid lines assume an opaque (f = 2.8) ionizing disk 
with (thin line) and without (thick line) a bremsstrahlung halo; 
the LMC contribution is shown by the long-dashed curve. The 
dot-dash curve is £ m (Ha) predicted for the upper side of the 
Stream due to the bremsstrahlung halo; the cosmic ionizing 
background is expected to dominate here. The open circles are 
the W 2 Ha measurements. 



The Stream Ha detections are (1.1, 0.63, 0.60) cm" 6 pc 
or, equivalently, (370, 210, 200) mR 2 . Fig. 3 reveals that the 
Stream is highly clumped: the non-detections and weak detec- 
tions fall in regions of low mean column density 1-3 x 10 18 
cm" 2 ), whereas the detections fall in regions of higher column 
density by at least an order of magnitude. The near-constancy 
of the detected signal argues for a photoionization origin. We 
suspect that the projected gas is clumpy on all scales, and that 
the detected regions are those in which the areal filling factor 
of neutral gas within the beam is essentially unity, so that all 
of the incident ionizing photons are absorbed by the gas. The 
close match between the Ha and Hi kinematics also supports 
static photoionization. The observed £ m thus directly measures 
the ionizing photon flux incident on this portion of the Stream. 
We can therefore use the Ha detections to determine the escape 
fraction of ionizing photons from the Galaxy. 

kpc (Feast & Walker 1987); for the SMC, (m-M) = 18.83 ±0.15, or 58.3±4.0 
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In Fig. 4, we show the predicted £ m values along the Stream, 
normalized to the Ha detections. This requires a mean opac- 
ity f ~ 2.8 ±0.4, corresponding to / esc ~ 6%. (Note that for 
f = 0, the Stream would be highly ionized.) The large opacity of 
the disk elongates the surfaces of constant photon flux density 
along the polar axis of the Galaxy (Bland-Hawthorn 1998). In 
this direction, the disk radiation field dominates over the cosmic 
background ($° < 4 x 10 4 phot cm" 2 s" 1 : Vogel et al. 1995) out 
to a distance of at least 200 kpc. We treat the LMC as a point 
source from which 5 x 10 51 ionizing photons escape per second 
(e.g., Smith, Cornett & Hill 1987). The Galactic radiation field 
dominates except within distances ~10 kpc of the LMC. 

4. DISCUSSION 

The conclusion that approximately 6% of the ionizing pho- 
tons produced in the Galaxy escapes entirely from the disk and 
halo has major implications for both the structure and ener- 
getics of the ISM and the cosmic ionizing background. The 
uncertainty in this estimate is approximately 50%, and is con- 
tributed equally by the uncertainty in the ionizing flux incident 
on the Stream (i.e., whether the beam filling factor of high col- 
umn density gas is actually unity for any of the measurements), 
and by our knowledge of the total ionizing photon luminosity 
of the Galaxy. 

Reynolds (1990) has estimated that approximately 15% of 
the ionizing photons produced in the Galactic disk are needed 
to power the diffuse ionized gas layer, i.e., they must escape 
from the disk and reach heights of several hundred parsecs be- 
fore absorption. Our inference that approximately 6% of the 



ionizing photon luminosity escape the halo entirely lends sub- 
stantial credence to the assertion of Miller & Cox (1993) that 
the O star population of the Galaxy is responsible for the ion- 
ization and heating of the diffuse ionized gas layer. 

A relatively intense, ionizing radiation field at large heights 
above the plane also has profound implications for the study of 
the high-velocity clouds (HVCs). Since their discovery more 
than 30 years ago, understanding of the nature and physical 
characteristics of the HVCs has been severely hampered by 
the almost complete absence of data on their distances (q.v., 
Wakker & van Woerden 1997). However, our determination 
of the ionizing radiation field in the Galactic halo, normalized 
by the Stream detections, indicates that it should be possible to 
derive Ha distances to HVCs out to at least ^100 kpc (Bland- 
Hawthorn et al. 1998; cf. Tufte, Reynolds & Haffner 1998). 

Finally, we note that if our value of / esc for the Milky Way 
is typical, then galaxies dominate over AGN in producing the 
cosmic ionizing background at low redshift (Giallongo et al. 
1997). 
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